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Abstract 

The performance of command and telemetry systems, useful in deep-space 
communications, is frequently affected by the radio-frequency phase error which is 
introduced at the point of reception by means of the carrier tracking loop. In 
low data rate communications, this phase error which is highly unpredictable may 
vary rapidly over the duration of the signaling interval. In this paper such 
phase variation is characterized by a sinusoidal input phase, k sin(u>ot +TT/6), 
which models a typical phase variation in communication over turbulent media. 
Conditions for synchronization stability and the acquisition behavior are examined 
by detailed computer study of the phase-plane trajectories for the second and 
third-order loops with perfect integrator. For k = 0.001, it is determined that 
f 0 /AK <_ 1/1» for system stability. Here AK is the loop gain. For given f 0 , the 
condition for stability is kf D <_ U, except for very small f Q . Thus the loop is 
jvlll useful under most fading conditions in deep space missions. 

It is demonstrated that for the phase variation considered the third-order 
loop has no real advantage over the second-order loop. Finally, it is shown that 
nonzero initial conditions may result in large steady-state phase error. 



Nonlinear Analysis of Phase-Locked Loops With 
Rapidly Varying Phase 
C.H. Chen and M. Fan 

I. Introduction 

The performance of command and telemetry systems, viseful in deep-space 
communications , is frequently affected by the radio-frequency phase error which 
is introduced at the point of reception by means of the carrier tracking loop. 

In low data rate communications, this phase error may vary rapidly over the 
duration of the signalling interval. Causes of this type of behavior in planetary 
entry are turbulence, dispersion, attenuation and residual doppler. The phase 
variations cannot be tracked by c. phase-locked loop of lower bandwidth, while the 
signal-to-noise ratio in this minimum loop bandwidth is too low. 

When the ratio of the system data rate to carrier tracking loop bandwidth is 
less than one, the problem of power allocation between the carrier and the data 
has been considered by Hayes and Lindsey tl] , Thomas [2J, Sergo and Hayes [3]. 
For channels with time-varying phase, Heller (U] examined the performance of a 
sequential decoding system. An excellent treatment of the nonlinear analysis of 
the phase-locked loops is given by Viterbi [5] and Lindsey [6] . 

In this paper the phase variation is characterized by a sinusoidal input 
phase, k sin(u Q t + ir/6), which models a typical phase variation in communication 
over turbulent media. Nonlinear analysis of the loop in the absence of noise is 
performed by extensive computer study of the phase-plane trajectories. Conditions 
for synchronization stability and the acquisition behvaior can be examined from 
the phase -plane analysis. Both the second-order loop and the third-order loop 
with perfect integrators are considered with zero as well as nonzero initial 
conditions. Comparison is also made between the second-and the third-order loops 
with rapidly varying phase. Some preliminary computer results were reported by 
Chen [ 7 ] and F an [8]. 
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II. The Loop Equations 

Following the notations of Viterbi ([ 5 ], Chapter 3), we consider first the 


differential equation of a second-order loop with perfect integrator. 


d 2 d 2 0 

— ® + AK cos* + aAK sin* = 

" Qw 


(1) 


dt dt 

where *(t) is the phase error, AK is the loop gain, 6^(t) is the phase of the 

input signal, and the transfer function of the loop filter is 
a 


F(s) = 1+2. 

s 


( 2 ) 


Ihe loop can track the frequency ramp with zero steady state error. Now we con- 
sider the important case that 0^(t) is varying several cycles over a bit interval 
of, say, 1 second which is typical in low data rate communications . The variation 


is normally caused by the time-varying channel. Let 
e i (t) = k sin(u^t + ) 

By normalizing the variables with 

2 

a' — - T ** AK+. 4* a ^ a" — ^ 

a A k * T AKt » ♦ dt * * .2 

d x 

Eq. (l) becomes g 

„ ku u x 

♦ + * cos* + a' sin* sin + £•) 

(AK) 2 AK 6 

which in turn can be written in the state equation form as 

2 


(3) 


( 4 ) 


(5) 


X 1 =X 2 


ku 

f — * - — 0 


0) T 


x 0 = - Xg cos - a f sin - 


(AK) 


• / O IT \ 

2 Sln hue * z > 


( 6 ) 


where x^ = *(t). It is noted from Eq. (5) that the larger the loop gain, or the 
loop bandwidth, the smaller the frequency of the forcing function given by Eq. (3). 
The frequency f Q is reduced by a factor of AK, and the amplitude ka) 2 is reduced 
by (AK) . In other words, the large loop gain reduces the effect of the time 
varying input phase *-(t). 

For a third order loop with loop-filter transfer function 
F(s) 


* a b 
s 


( 7 ) 
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the differential equation (Viterbi [51, p. 64) is 

d 2 f t a 2 g 

^4 + (AK + aAK)sin4»(t) + bAK sin*(u)du = =■ (8) 

, , c at; J 

dt ' dt 

o 

which, using Eqs. {k) and (3), can be reduced to^ 

„ t kw (i) IT 

$ + $> cos$ + a* sin$ + b* f sin<fr dx = - — — 0 sin (ttt- + £) (9) 

i (AK ) 2 M 6 


which in the state equation form becomes 


*1 = X 2 


cos x^ - a f sm x 


1 - b- | si 


ko»: 


U> T 

O . TT, 


sin x, dt Q sin (— jrp + ■?) (10) 

(ak) 2 AK 6 


where 


b' = and x = $(t). 

(AK) 2 1 


Phase-plane analysis of Eqs . (6) and (10 ) is performed by using the second- 
order Runge-Kutta method (see, e.g. [ 91 ). The computer results are reported in 
the following sections . 

HI. Nonlinear Analysis of the Second-Order Loop 

Consider first k s 0.001. The loop behavior depends on the ratio f Q /AK. 

Let a* * 1/2 and AK = 32, the phase-plane plots of the loop are shown in Figs. 1(a), 
(b) and (c) with f = 6.4, 8.0 and 16.0 respectively. f Q = 8.0 appears to be the 
threshold value above which the loop would not be able to track the input phase. 

For f Q /AK 1/4 the loop would settle with a stable "limit cycle." It is noted 
that the steady state error cannot be reduced to zero because of the continuous 
input phase variation. Thu3 f Q /Ak < 1/4 is the condition for stability. By 
increasing the nrnber of points to 3000, Fig. 1(d) shows that the steady state 
trajectory drifts only slowly. The loop reaches the stable trajectory after less 
than two cycles of change in phase. 

Next we examine the threshold value of k for specified f Q . Figs. 2, 3, 4, 5, 6 
are the phase-plane plots for f = 1.1, 1.0, 0.8, 0.5, 0.25 respectively. The 



threshold here is not a critical value but it may be concluded that kf s£ 4 is 

o 

the required condition for stability with AK a l6, a = 8. When f becomes small 
the allowable product kf Q rises rapidly. In the limit f Q - 0, kf Q can be any 
positive value. If the phase error varies one full cycle in one second, the 
above condition states that the maximum allowable phase error is 4 radians . This 
indicates that the loop can function properly under most fading conditions in 
deep space missions. 

For example, based on the updated knowledge of the effect of turbulence in 
the Venus atmosphere on radio propagation [101, we anticipate that the phase-locked 
loop with proper bandwidth can maintain a continuous communication between the 
atmospheric probe and the Earth. 

IV. Nonlinear Analysis of the Third-Order Loop 

Consider the case k = 0.001, AK = 16, a = 8, b' = 1/32. Figs. 7(a), (b), (c) 
are the phase-plane plots of $ vs $ for f = 3.2, 4.0 and ' .0 respectively. For 
the same ratios of f Q /AK, the trajectories of the third-order loop drift faster. 
f Q /AK <_ 1/4 also appears to be the threshold condition for stability. The increase 
of b 1 to b* = 1/16 as shown in Figs. 8(a), (b), (c) only causes the trajectories 
to drift more. The effect of increasing k is clearly illustrated in Figs. 9(a)-(d). 
The threshold value for k is in the region 0.1 < k < 1.0, The condition for 
stability is kf o < 4 which is consistent with the second-order loop results. A 
careful comparison is made between the second-order and the third-order loops by 
using exactly the same loop parameters AK = 8, a = 4, k = 0.001, f = 1.0 with 
printout of 1000 points. Both loops are stable and take the same amount of time 
(50 seconds) to cover 6 cycles of trajectorits. The results are shown in Fig. 10(a) 
for the second-order loop and Figs. 10(b) and (c) for the third-order loop with 
b* * 1/32 and 1/8 respectively. It is noted that the third-order loop has larger 
phase error which increases with b' . Based on the above study, we may conclude that 
for the sinusoidal phase variation considered in this paper, the third-order loop has 
no real advantage over the second-order loop. 
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V. Effects of Nonzero Initial Conditions 

To study the loop behaviors under different initial conditions, we consider 
both the second-order and third-order loops with AK = 8, a = 4, k = 0.001 for 
the initial conditions ($, «>) = (-3.14, 6), (-3.14, 4), (-3.14, 2), (-3.14, 0). 
The results are shown in Fig. 11 for the second-order and Fig. 12 for the third- 
order loops with b' = 1/32. The trajectories in the lower half plane are the 
mirror image of those in the upper half plane with respect to the <j> = 0 axis. 

It is interesting to note that, in all initial conditions considered, the loops 
reach a steady-state of $ = 0 and $ = constant in spite of the sinusoidal input 
phase. Hie steady-state phase errors may be too large, however, from practical 
viewpoint. Thus the nonzero initial conditions should be avoided if the large 
steady -state phase error is not tolerable. 

VI. Concluding Remarks 

We have examined the critical parameter values of the second and the third- 
order loops with rapidly varying phase which is modelled by a sinusoidal input 
phase variation. The nonlinear analysis is performed by studying the phase-plane 
trajectories. Although the third-order loop does not have the advantages as it 
normally has in tracking phase errors, some modification of the loop structure 
appears necessary to obtain a more efficient tracking system. Methods of signal 
acquisition aids, as suggested by Lindsey [6], especially the computer-aided 
acquisition for low signal-to-noise ratios should also be examined. 
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